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Abstract 



We apply the concepts of multifractal physics to financial time series in 
order to characterize the onset of crash for the Standard & Poor's 500 stock 
index x{t). It is found that within the framework of multifractality, the "anal- 
ogous" specific heat of the S&P500 discrete price index displays a shoulder 
to the right of the main peak for low values of time lags. On decreasing 
T, the presence of the shoulder is a consequence of the peaked, temporal 
x{t + T) — x{t) fluctuations in this regime. For large time lags (T > 80), we 
have found that Cq displays typical features of a classical phase transition at 
a critical point. An example of such dynamic phase transition in a simple 
economic model system, based on a mapping with multifractality phenom- 
ena in random multiplicative processes, is also presented by applying former 
results obtained with a continuous probability theory for describing scaling 
measures. 
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I. INTRODUCTION 



In this work we apply the concepts of multifractal physics to financial time series in 
order to characterize the onset of crash for the Standard & Poor's 500 stock index. We shall 
present an example of dynamic phase transition in a simple economic model system based 
on a mapping with multifractality phenomena in random multiplicative processes and by 
applying former results obtained with a continuous probability theory for describing scaling 
measures. 

An attempt is made to characterize the presence of stock market crashes by solving for 
a price equation from a non-linear equilibrium model and showing how multifractal physics 
measures can be generated from this equation. We found that an " analogous" specific heat 
Cq of the S&P500 price data displays a shoulder to the right of the main peak as a function 
of time lags. For large time lags, Cq resembles a classical phase transition at a critical point. 
We explain this dynamic phase transition by a mapping with multifractality phenomena in 
random multiplicative processes. Within this description the temporal price variations of a 
commodity displays features of an "analogous" phase transition from infiated to devalued 
prices, when the excess demand is not linear in the asset price. An analytical expression for 
Cq of the economic model system is derived. 

Finance and physics joint "ventures" have attracted considerable interest in the literature 
for over many years [|1|. These efforts have allowed to pursue analogies between stock 
market dynamics and stochastic models commonly used in the statistical physics of complex 
systems . Such parallel analysis have been useful to best quantify and understand possible 
correlations in financial data by measuring the autocorrelation function and the power 
spectrum 0]. Another notable example is the analogy with the scaling properties observed 
in turbulence 0-0]. 

Parallels have also been drawn between a very simple theory of financial markets and 
the quantum gauge theory ^ (see also comments in p). Other interesting studies for 
investment strategies in diversified portfolio stocks have been recently discussed in terms of 
products of random matrices and of multiplicative random walks [jll|. Realistic price 
fiuctuations have also been found to emerge in the adaptation of a system {i.e., the traders) 
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to complex environments |jT2|] and from the self-organization of a (closed) system of traders 
without external influences p3|] . 

Thus, the financial market dynamics is still an open subject for physicists and economists. 
In particular, the basic principles governing the origin of stock market crashes are far from 
well understood by both communities , specially in regard to the problem of world-wide 



market dynamics. Using the renormalization group theory []T^|-|T^ , it has been proposed that 
this cooperative phenomenon could be the result of a critical phase transition (see also [|T8[] ). 
Motivated by these suggestions, it is then tempting to determine the constrains needed to 
understand and describe analytically an analogous phase transition from a new perspective. 
That is, by using the characterization of multifractal singularities where an analogy with 



thermodynamics has already been established [0,0 . 

Multifractality was initially proposed to treat turbulence and, in recent years, it has been 
applied successfully in many different fields ranging from model systems such as Diffusion 
Limited Aggregation to physiological data such as heartbeat. The multifractal analysis of 
financial discrete sequences developed in our paper is another aspect in a relatively new 
topic of physics, the so-called econophysics. 

II. MULTIFRACTAL CHARACTERIZATION 

This paper will focus on scaling laws in financial returns, particularly on the absolute 



values of returns and their scatter as a function of time. Similarly to other studies [0], we 
are not attempting to predict a price drop or rise on a specific time on the basis of past 
records, but we are after a new characterization of the presence of stock crashes, which is 
the added value of our distinctive approach. This is important to mention since forecasting 
has been extensively studied in the econometrics literature by other techniques such as 
ARCH and multivariance models depending, e.g., on the data seasonal behavior and previous 
experiences of the forecasters. Interesting work on fractional processes of these types in 
econometrics is very well known and these are beyond the scope of our work. Important to 
mention is that recent theoretical work has considered the possibility of fractal nature of the 



absolute value of returns on exchange rates |21 



Let us introduce next our new characterization of a financial time series x{t) based on 



the well-known definitions used in multifractal physics. In Fig.l we show the temporal 



x{t + T) — x{t) fluctuations of the S&P500 index (data available in |^) for the 12 years 

period 2 January 1980 - 31 December 1992 as a function of the trading time lags 1 < T < 220. 

From this figure, it can be observed that the peaked narrow fluctuations measured for T = 1 

spread out on increasing T. For T 80, the maximum and minimun difference values of 

the so-called "Black Monday" crash measured in 1987 (shown as the largest straight line for 

T = 1 in the figure), become comparable to the valleys and picks differences of the relative 

S&P500 fluctuations over wide periods of time. This peculiar behavior is shown next to lead 

to an analogous thermodynamic phase transition when varying T. 

Let us consider the following measure over N intervals 

\x{t + T)-x{t)\ 
'~El,\x{t + T)-x{t)\ ' 
with T representing some finite trading time lags T. Clearly, the above relation can be 

viewed as a normalized probability measure with /ij > 0. 

From fii in Eq.(|^), we then construct the corresponding generating function Z, and its 

moments q, which follows the scaling 

N 

Ziq,N)=\unY.^^^iTr-N-<'^ ■ (2) 



t=i 



To get a thermodynamic interpretation of multifractality (see, e.g., p!9|j20|j25[] ) , we divide 
the one-dimensional system of length L into lines of length i; thus ~ L/i. We then 
associated this with the number of discrete x{t) time sequences considered in Eq.(0) in 
order to relate T, L and i in the definition of the measure. For i/L ^ 0, the function r 
relates to the generalized fractal dimensions {q — l)Dq. Similar multifractal analysis has 
recently been performed for the energy dissipation field of turbulence logistic maps ||24 
and surface roughening p5| . 



By following the thermodynamic formulation of multifractal measures, we can also derive 
an expression for the "analogous" specific heat as follows 

C,^-^^^r(g + l)-2r(g)+r(g-l) . (3) 

For large T, we shall show that the form of Cg resembles a classical phase transition at a 
critical point. 

We shall return to these equations in Section III. 



III. DYNAMICAL MODEL 



Similarly to ||T3| , p6| in our approach there is only one stock. To study the price changes 
for one commodity it is necessary to derive a dynamical equation which results from the pre- 
vailing market conditions. The market is usually considered competitive so it self-organizes 
to determine the behaviour of prices. We assume here that all factors determining the 
demand D and the supply Q other than the asset price p remain constant over time and 
denote these quantities in equilibrium with an asterisks (*). In the following all variables 
are dimensionless. 

In a competitive market it is expected that the rate of price increase should be a functional 
of the excess demand function E{p) = D{j>) — Qijp). Hence one writes dp/dt = f[E{j))] p7| . 
Assuming that a commodity can be stored then, in general, the flow of demand does not 
equal the flow of Q output. Hence stocks of the commodity (or product) build up when the 
flow of output exceeds the flow of demand and vice-versa. Then the rate at which the level 
of stocks S changes can be approximated as dS/dt = Q{p) — D{p). From these relations, 
a price adjustment relation that takes into account deviations of the stock level S above 
certain optimal level So (to meet any demand reasonably quickly) is then given by 

+ . (4) 

where 7 {i.e., the inverse of excess demand required to move prices one unity pGf) and A 
are positive parameters. If A = 0, the price adjusts at a rate proportional to the rate at 
which stocks are either raising or running down. If A > 0, prices would increase when stock 
levels are low and rise when they are high (with respect to So). We shall assume here A to 
characterize a noise term. 

In our description, for each asset price p, we postulate simple non-linear forms for the 
quantities D demanded and Q supplied such that 

D{p) =d*+do[i- ^^{p-p*r +...]{p-p*) , 

Q{p)=q* + qo[l-^^{p-p*? + ...]{p-p*) , (5) 

where do, qo and d* = D{p*), q* = Q{p*) are arbitrary coefficients (related to material costs, 
wage rate, etc), p* is an equilibrium price and 5 is our order parameter. We write D and 
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Q as a Taylor series expansion with the usual linear dependence (independent of S) plus a 
non-linear correction. Higher order terms 0{A) are here neglected for small p — p*. In the 
above we might also consider two different 6s, but to reduce variables to a minimum we 
assume D and Q to vary similarly from linearity. To simplify notation we also define 

Po = qo-do . (6) 



In the context of a simple economic model , it is reasonable to assume that So depends 
linearly on the demand; e.g., So = io + ^-D, with io a constant and i satisfying the constrain 
below. The postulated linear dependence of the optimal stock level So on D at equilibrium 
provides a complete economic model as in ||2^. Therefore, in equilibrium (where ^\p* = 
and ^1^. = 0, so that demand equals supply and S = So), from the above we obtain 

d* -q* = , S* = io + i{d* + doP*) . (7) 

And after some little algebra we find that the price of one strategic commodity is governed 
by the general equation 

A,3.(P-P*)|1- %{p-f?]^0 . (8) 



The linear case is for 5 = 0. This leads to p{t) — p* (x Ai cos{t^J\j3o) + A2sin{t^J\l3o). 

A. The Simple Case IXdo = 'yPo 
To keep the mathematics simple we choose p* = and consider i to satisfy 

\do 

The general Eq.(|^) then reduces to 



(9) 



^ + A/?oP-^^p'~0 . (10) 

This is our dimensionless price adjustment equation which gives rise to a burst as discussed 
later. When 5 7^ and [XPo, ^^A/3o/2] > 0, it has the well-known kink solutions 
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Pit) = ±^ tanhi^ t) . (11) 

Clearly jSo of Eq.(^) must be positive. Since in a free market economy the demand for a 
product (or commodity) fall when its price increases, then it is reasonable to assume do < 
in Eq.(ll). As the price raises, the supply also increases; hence in general one also assumes 
qo > 0. These conditions yield /3o > as requested in Eq . (|TTD and also doi > 0. 

In the case 6 < the above function for p{t) displays a sudden decline around the 
equilibrium value p* taken to be pit*) = p{0) = 0. 

B. The Case IXdo + -ijio 

If we consider the case in which iXd^ ^ 7/5o, the price equation of Eq.(|) results in a 
Lienard-type of equation: p" + gi{p)p' + go{p) = due to the presence of the ^-term. What 
we present next is a brief discussion regarding its possible solutions. Again, we set p* = 0. 

With the aid of the substitution ^ = w{p), so that ^ = w{p)^, the Lienard equation 
can be reducible to an Abel equation of the second kind ww' = fi{p)w + f^ijp)- For small 5, 
the substitution w{p) = p^K{p) + ^{iXdo — 'jPo) leads to the Bernoulli equation with respect 
to p = p{K): 

3Kip)[^jiiXdo - 7/jo) - Kip)]^ = K{p)p + ^^"^ , (12) 

whose solution is 

= [jiiXdo - 7/5o) - Kip)]-'/' { Ko + X 

""^mp)^ ^^^^3'3'3'¥ K{p) ' 
with Ko an arbitrary constant and 2-F1 the first, complex hypergeometric function as arising 
in many physical problems. It converges within the unit circle \{iXdo — 'jjSo)/ K{p)\ < 2/5^. 

By examining these equations one recognizes that if iXdo 7^ 7/3o, other behaviour might 
appear for p{t) (different from the one in Eq. (pJ]) ) . However, it is important to note that 
such possible behaviour can essentially be found in the complex plane since the solutions of 
Eq.(p!3D are driven by the i^~^/^-term and depend whether 7/^0 is greater or smaller than 
iXdo- A detailed analysis of such solutions is beyond the scope of this work. We only study 



here the simplest dynamical economic model where all its variables, including prices, are 
related to the demand and supply functions as seen in a free market. 



IV. MODEL AND ANALOGOUS PHASE TRANSITION 

Let us now identify the behaviour of p{t) with an analogous phase transition as seen 
in multifractals. Following the analogy with critical phenomena in the time domain as 



proposed in []T6| , we also consider t to be the relevant variable for the analysis of an possible 
existence of an analogous critical point. To derive a connection between our economic model 
and multifractality we first briefly review multifractal phenomena in Random Multiplicative 
Processes. 

A. Multifractality in Random Multiplicative Processes 

Multifractality emerges in random multiplicative processes for a self-similar function i/j 
such that is rescaled as ip^x) = e~^^^^^ip{x), where L{1) is the generalized Lyapunov exponent 
for the first moment of ip and x is a space variable for A^-disorder fluctuations on a unit 



interval |^. All information about these systems is embodied in the non-zero, positive %jj 
measures. 

To analyse the analogy of multifractality with thermodynamics one then scales the mo- 
ments q of the functions ip with respect to segments / < 1 as 

Z(g,/)= hm E^^,^~r(^) , (14) 

fc=l 

which defines the exponents r and Z is a formal partition function. 



It has been previously shown that for general random multiplicative processes |28|, r 
satisfies 

r(g) = -(l-g)- - gL(l)} , (15) 

where 

1 ' 
fc=i 

h= lim — In/-^ . (16) 



We shall use these findings to derive a connection between our economic model variables 
and multifractality phenomena. From this connection, we shall identify all the model quan- 
tities that might exhibit multifractality within the framework of a stochastic multiplicative 
process. These processes are known to generate power law probability density functions 

2|,|3g. 



B. General Continuous Probability Theory 



We use next the simple continuous probability theory discussed in that allows to 
explore the genesis of an "analogous" phase transition from the point of view of a nonlinear 
singularity spectrum equivalent to multifractals. A crucial feature of this formalism is to 



consider t\J A/3o/2 to be a continuous random variable. Then within the framework of general 
probability theory (see, e.g.^ [Hi); the distribution function of this random variable, defined 
in a line and in terms of its probability distribution P, can be approximated as P{Ci < C' ^ 
C2} = ^((2) - ^(Ci) ~ /c? <PiC) dC, where {} indicates the funct ion interval and 6 is a 
uniform probability density that needs to be specified. 

Inspired by well-known results for the probability distribution function of real economic 
data with power-law tails p9| , |30| , |33| -p5[] , we assume 0(C) = (0o/2)[l — (5 /\^)H{Q] such 



that 0(C — > +oc) — > and —00) ^ 0o > 0, with H given by the real solutions of 

a static, dimensionless Ginzburg-Landau-like equation. It resembles the spin-fiip function 
of the well-known Glauber-Ising chain model. Using such solutions, we shall show next 
that it is possible to establish a relation with thermodynamics similarly to multifractality 
phenomena. 

It is at this point that we relate H{() to the positive solutions p{t) of Eq.(|TUp and map 
C/Co^ ty'A/So/S within the framework of general probability theory. Following Ref. |]HT]], it 
is then straightforward to derive an expression for an "analogous" specific heat, , for our 
economic system. To obtain such an expression we evaluate first the integral of P over the 
range [Co,C]- Thus we have 

g{0-Q{Co) = ^j\l-^tanh{^J}dC^T{0 , (17) 
which, in turn, defines the function t{(). The Q functions satisfy Q{Q > g{(o), or alterna- 



tively, 0(C) > 0(Co) (since dQ jdC, = 0(C) ^^), which is in accord with the above assumption 
for 0(C ±00). 

Let us see next how our analysis from general probability theory would capture multi- 
fractality. Similarly to the dielectric breakdown model or the Poisson growth model, where 
the local field is set proportional to the growth probability |3^, we assume here that the 
continuous function r (to be identified as a free energy) is proportional to the probability 
distribution P as in Eq.(|l7|). 

After a little algebra, the above integral gives 

r(C/Co) ^ (1 - C/Co)r(0) - ^{Inco.MC/Co) - (C/Co) Incos^l)} , (18) 

in which r(0) = ^(0) - GiQ = -^{1 + -j^^x}, such that Tx = - lncosh{l). We have 
assumed that 5Tx/^/2 << 1, hence 6 < 0. 

Up to this point we have not carried out any actual numeration or scaling of a particular 
fractal configuration or set as used in multifractal theory. However, to gain further insight 
and explain how the present general probability can mimics multifractal behaviour, let us 
now relate the above r to a measure within a random multiplicative process. By mapping 
our results to such a measure, all other quantities that can be scaling and are directly derived 
from r, such as an analogous specific heat, will then follow. 

From a comparison between Eqs.(^) and (pisj) one can easily identify the following terms 

h ^ 2V2 ' 
t(0)^-1 , 

L(C/Co)^lncosh(C/Co) . (19) 

It is from this simple mapping between our results and those for random multiplicative pro- 
cesses that we can made an attempt to understand how the existance of the economic model 
order-parameter 6 (and from it, the non-linearities in the demand and supply functions) 
leads to obtain multifractal-like behaviour. 

From the above mapping we deduce that if 5 ^ 0, then the moments q of the functions 
ip with respect to segments / would vanish since / ^ 0. It follows that a linear economic 
model would never exhibit multifractal features since in this case r(g) = q — 1. Furthermore, 
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it is important to note that a Lyapunov exponent of the type we derive resembles that of 
a random multiphcative process with A described by the probabihty distribution -P(A) = 
(1/n) YJi=i — ^i) By considering the actual definition of L{q), the moments of 

the exponential measures in a random multiplicative process might well be related to our 
reduced variable as g ^ C/Co ^ t^A/9o/2. 

C. Possible Multifractal Features 

To analyse multifractal features in our economic model, when mapped into a random 
multiplicative process as discussed above, we consider standard definitions: t{(/(o) = 



[(C/Co) — 1]-Dc- According to such definitions [p^-|iO[], the function r represents an "analo- 
gous" free energy and D^^ the multifractal dimension. 
From Eq. ([T8|) it follows that 

Dc ^ -r{0) + ^7f^^^^{l^^«^MC/Co) - (C/Co) \ncosh{l)} , (20) 

such that ( 7^ Co- From this relation we obtain -D^^o = ~t(0), and -D(^^+oo = -D^^o ~ 
+ r,}; D^^_^ = Dc_o + - r,}. If C = Co, then D^^^^ = D^^o - %f {Fa + 

tanh{l)}. 

Complementary to r we also define 

"(C/Co) = d(^f(0 ^ D^^o - ^{Fa + tanh{C/Q} . (21) 
It can be easily shown that a^ax = a(C/Co ^ -oo) = D^/(^^^_^, and amin = «(C/Co ^ 

+ 00) = L'^/f^^+oo- 

Also according to multifractality phenomena, a possible analogy with thermodynamics 
can be established by relating r to /(a) = (C/Co)tt(C/Co) — '^{(/O via a Legendre transfor- 
mation. From this analogy, where / is athe analogous "entropy", our analytical expression 
for the analogous "specific heat" of the economic system becames 
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V. DISCUSSION 



Let us see next how the type of behaviour given in Fig.l influences the "analogous" 
specific heat function Cq deflned in Eq.(^ and how it might characterize the onset of a 
crash for a real stock market index. In Fig.2 we plot the "analogous" specific heat Cq of 
the S&P500 index for four different time lags T = 1, 10, 30, 80. The + + + curve is for the 
1984-1988 data, xxx is for 1982-1990 and ooo is for 1980-1992. 

For time lags T — > 80, we find that the main peak of our numerical Cq resembles a 
classical (first-order) physics phase transition at a critical point given by the main peak 
position. The peak turns symmetric around the value g = — 1. Surprisely, this "analogous" 
specific heat Cq of the S&P500 index also displays a shoulder to the right of the main peak 
as a function of smaller time lags. Clearly, on decreasing T, the presence of the shoulder is 
a consequence of the large, temporal x{t + T) — x{t) fiuctuations in this regime. We note 
that such peculiar behaviour for a double peaked specific heat function is known to appear 



in the Hubbard model within the weak-to- strong coupling regime pi |. 

The relation in Eq.(0) requires T — ^ oo where the shoulder tends to vanish. It is this 
feature that make us believe that a large crash for an stock market index can be characterized 
by an " analogous" specific heat which resembles a classical phase transition at a critical point 
as studied in multifractal physics. 

We now turn to the results of our theoretical economic approach. For T = 80, the full 
line in Fig.2 represents theoretical results for from Eq.(P^) by choosing ^ g -|- 1 to fit 
the main peak position. It can be seen that the theoretical curve resembles the phase 
transition features of the real S&P500 economic data. Our simple model for one commodity 
shade light into the main observed features regarding a possible analogous phase transition 
that occurs when the excess demand becomes non-linear (c./., cubic p-term in Eq. (p^ ) in 
terms of the price for one commodity. We believe the width difference between both curves, 
i.e. the analytical and the estimated C{q) function from historical S&P500 data, is due 
to the fact that the S&P500 price index is made of large-capitalization stocks representing 
a "basket" or portfolio of commodities. 

For large time lags, there is a sharp peak that resembles the quantitative signals measured 
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in multifractals p8|-^. From this feature we presume the existence of an analogous, say, 
critical point above which inflated prices for one strategic commodity might be found. 
The maximum and minimum values of a in Eq.(PTD (for more details see also |^) allows for 
the existence of a critical point (* above which the infinite hierarchy of phases can be found, 
but below which a single phase appears characterized by amax- It resembles a classical phase 
transition at a critical point \TE\. 

Of course, the analogy between multifractality and a thermodynamic phase transition 
as discussed here does not imply that the economic system has a phase transition. What 
we have shown, as a direct consequence of the p^-term in Eq . ([T0|) , is that prices can became 
inflated prior to equilibrium (i.e., t < by convention), whereas after a sudden crash prices 
might devalue. The greater 1 5 1 -values are taken, the smaller the price reduction becomes. 
If p > the amount supplied exceeds demand and stocks accumulate whereas if p < the 
stocks deplete. Also by tuning 5 — in Eq. ([TT|) , we are able to predict that prices can 
decrease (or increase if > 0) monotonically with C — > for all (. 

It is a well-known fact the hill behaviour of the generalized dimensions Dq for g < 



when using the box-counting method as in the present work p2| , |43|] . Thus, we check if Dg is 
sufficiently smooth for Cq to be meaningful. In Fig. 3 we represent Dq of the S&P500 index 
for the 1980-1992 period within the range q > for the time lags T = 1, 10, 30, 80, 120. The 
total number of sequence points analysed include 1495 points for the different Cq curves, 
1500 points for Dq and 9864 points for the S&P500 1980-1992 data set. We also fit this 
non-seasonal data by standard exponential smoothing techniques, and from the fitting our 
error estimates for the Dq curves is found to be less than 3%. Our results for Dq at negative 



q (not shown) follow a typical convergent behaviour as can be seen, for example, in |31|. 

From Fig. 3 it can be seen that, when increasing T > 30, Dq is fully mult ifractal- like and 
for T > 120 it becomes fiatter {i.e., uniform measure) and tends to one independently of q, 
so multifractality becomes smaller. For lower values of T < 30 we find a non-monotonous 
decreasing behaviour of Dq, conceivable within the double peaked form of Cq displayed in 
Fig. 2, which relates to the presence of the onset of crash for the S&P500 stock index in Fig.l. 
Using this data set, we also estimated the multifractality strength of the time sequence by 
considering the limit 1 — Dq^^o for the few different time lags displayed in Fig.3. We find 
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that that this quantity does not follow a power-law scaling for T as in citeRom97 for values 
T ^ 1. This is a consequence of the complex network of trading interactions comprised in 
our non-linear forms for the supply and demand functions. 

The present choice for a excess demand function of the form E{p) oc {p ~ p*) — ^{p — p*)^ , 
with 5 7^ plays a key role in our description. We have d* = q*, hence the second price 
derivative ~ SPo^^P (or that of D and Q) is price dependent. It is such a behaviour of 
E (independently of the sign of 6), that leads to obtain an abrupt fluctuation in the price 
dynamics and characterize multifractality phenomena. 

Our expressions for the demand and supply functions of Eq.(^ are justified as follows. 
As seen in Fig.4, the (commonly used) linear p-dependence for D and Q and our assumed 
non- linear form for the these functions display similar behaviour when \6p\ << 1. Even 
more important, this figure depicts the fact that as price falls, the quantity demanded for 
a commodity can increase in agreement with one of the basic principles of economy. Since 
the demand curve can indeed change in a number of ways which may not be at all obvious. 



similar tails to our D and Q functions has also been previously hypothesized in ||44|| . 

In real world, exceptions to the general law of demand can take place making the D curve 
to increase upwards from low- to high-prices. However these exceptions, which include goods 



of conspicuous consumption -as, e.g., certain articles of jewelry- are not very important ps 



Theoretically, this simply would mean to set do > Oin our D{p) function independently of the 
order parameter 6. The upward dependence of D can occur as soon as there is speculation 
(as in precious stones mentioned but also for market prices). The idea is that when price 
increases, investors buy, because they hope the price will keep climbing. This is called a 
"trend-following" investment strategy. On the other hand, our choice for Q (with qo > 0) 
also follows the typical behaviour observed in a competitive market (where no individual 
producer can set his own desired price). That is, the higher the price, the higher the profit, 
then the higher the supply. 

It is also important to mention that the present choice for D and Q, and their linear 
p-dependence both lead to the same linear relation for D vs. Q, namely {D — d*)/do = 
{Q — q*)/qo- Such a linear behaviour founds frequent use in applied economics. 

Furthermore, it is well known in econometrics that demand functions are somewhat 
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abstract quantities since all of these data are taken to refer to possible events at just one 
moment of time [^. In particular, consumer behaviour [i.e. tastes, desires, . . .) can shift 
a (concave up or concave down) D curve (which may resemble the tails of our D and Q 
functions). Also, typical demand and supply of capital look like a step function underlying 
forces toward interest rates |^ . 



Usually additional hypothetical information is needed to make up realistic demand re- 
lationships {e.g., consumer interviews). In principle, one might also determinate D out of 
many data sources from different economic sectors using standard statistical techniques for 
multiple regression time series analysis. But the effectiveness of such an approach varies case 
by case. In view of all of this, our non-linear approach for D and Q may well be placed for 
simulating market situations leading to a sudden decline around equilibrium for the price of 
a commodity. To this end, we add that a commodity price function displaying an inflection 
point at a characteristic frequency has also been theoretically discussed in [^]. 

Another important test for our D and Q expressions arises by considering possible ag- 
gregated changes in conditions of demand (or, supply). Because of such aggregated changes, 
e.g. due to buyers' income and scale of preferences, it is always difficult to estimate how 
much of the change in D is due to price alone. To know what these effects are upon D, and 
obtain the degree of responsiveness of demand to price variations, it is necessary to study the 
point price elasticity of demand curves, defined as e = —dlogD/dlogp This quantity 



is said to be elastic (or fiat) if £ > 1 or inelastic if £ < 1 depending on the demand schedule. 

By assuming a linearly decreasing p-dependence for D one obtains the expression 
^ = - i+idhdoP) - if 5 ^ one gets e ^ - [..^^slrJ^i+ZJ /^^^y Hence we find that a 

transition from an inelastic to an elastic demand curve appears at p ^ — i_^Jp/°2)^ Whereas 
by considering the common linear dependence of D{p) one finds p ~ —d* /do. Therefore, 
both approaches lead similar results for the elasticity of demand (or supply) if 5 ^ 0. 

VI. CONCLUDING REMARKS 

We have found that within the framework of multifractal physics, the " analogous" specific 
heat of the S&P500 discrete price index displays a shoulder to the right of the main peak 
for low values of time lags. On decreasing T, the presence of the shoulder is a consequence 
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of the peaked, temporal x{t + T) — x{t) fluctuations in this regime. For large time lags 
(T > 80), we have found that Cg displays typical features of a classical phase transition at 
a critical point according to multifractal physics. 

Our simple continuos model for one commodity mimics the main observed features of 
C{q) for large trading time lags. We believe the width difference between the analytical 
and estimated C{q) curves from historical S&P500 data is due to the fact that the S&P500 
index comprises many commodities. From these results we conclude that an analogous phase 
transition might occurs when the excess demand becomes non-linear (c./., cubic p-term in 
Eq.(p!oD) in terms of the commodity price. 

We have assumed that there is only one stock in which a commodity can be stored. The 
market has been considered competitive so it self-organizes to determine the behaviour of 
prices. All factors determining D and Q other than p are assumed to remain constant over 
time. There exists a price adjustment relation that takes into account deviations of the 
stock level 5* above certain optimal level So characterized by a noisy A parameter. We have 
postulated simple non-linear forms for the quantities D demanded and Q supplied, with 5 
the order parameter, and have neglected higher order terms in their expansion on p — p*. We 
have followed the context of other simple economic models and assumed that the optimal 
stock level 5*0 depends linearly on the demand (with i being the slope). We have considered 
i to satisfy the constrain in Eq.(^, relating the economic model variables: 7, A, qo and do- 
We have d* = q* and set the equilibrium price p* to 0. We have identified the behaviour of 
p{t) with an analogous phase transition as seen in multifractals, by considering t to be the 
relevant variable and using a simple continuous probability theory. 

We have related t^\(3o/1 to a continuous random variable and have related its prob- 
ability distribution function to our solutions for p{t) given in Eq.(|l^). Our definitions for 
the analogous thermodynamic variables have been done according to definitions used in 
multifractal physics and by mapping to multifractality phenomena in random multiplicative 
processes. 

Of course the scenario of a transition from, say, inflated to devalued price changes in 
the time domain is pure speculation. However a great deal of relevant information has been 
extracted from the present continuous approach which, essentially, does relay on 5 only. 
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Our description presumes the existence of a stationary probability function P as in Eg. (p!7|) . 
We have assumed P to be proportional to the continuous function r (identified as a free 
energy) similarly to the dielectric breakdown model or the Poisson fractal growth models 
p6| . We add that this type of approximation is also used when modeling earthquakes, where 
the probability function of the total number of relaxations (size) of the earthquakes is set 
proportional to the energy release during an earthquake . 

Within the framework of general probability theory (see, e.g.^ HH); a continuous arbi- 
trary random variable can have an associated probability distribution (if it is discrete) or 
a positive probability density (if it is continuously distributed). The later is required to be 
differentiable. These can then be related by an integral equation of the type used in this 
work to then derive the analogous thermodynamics equations. To do this we have proceed 
as follows. We have used the characterization of multifractal singularities -where an analogy 
with thermodynamics has been established in the literature- to derive, and associate a mean- 
ing to our analogous quantities. The main point to understand is that we have not carried 
out any direct numeration or scaling of particular fractal configurations or sets, but we have 
evaluated the integral of the probability distribution P (for the continuous random variable 
Q assumed to be related with our p{t) function and mapped the results for r to those of a 
random multiplicative process. Our equivalent definitions for the thermodynamic variables, 
derived by solving such integral, followed from the convention used in multifractal physics in 



the sense that r must not be a linear function of C (see, e.g., [^)- Alternative approaches 
can also be found in l5^j5l| . In all theses cases, analogous thermodynamic quantities follow 
from the Legendre transform of t{Q. 

As discussed in [^,^], the concept of a phase transition in multifractal spectra was 
first found in the study of logistic maps, julia sets and other simple systems. Evidence was 
then found for a phase transition in more complex random systems such as diffusion limited 
aggregation. The condition for an analogy between multifractality and thermodynamic phase 
transition is that the analogous free energy (a in our notation) undergoes a quite sharp jump 
near a critical value C* . For values of C < C*5 the analogous free energy r is dominated by 
the maximum energy term Umax and a singular behaviour of the specific heat is found at 
this point. A well-defined analogous entropy function f(a) also suggests the existence of 
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an analogous phase transition. We have shown that the present kinetic description satisfies 
such properties obeying the constrains for £ in Eq.(^. This condition simphfies the analysis 
and corresponds to the case in which the term for the first derivative of p with respect to t 
{i.e., ^) is absent in the price adjustment Eg. (]T0|) . 

Application of multifractal analysis to discrete ID time sequences as derived, for example, 
from the cellular automaton model of a rice model is no new (see, e.g., |^). Nevertheless, 
to our best knowledge, we have related multifractal physics to financial time series for the 
first time. Our main contribution has been the analysis of the "analogous" specific heat (or 
second derivatives) Cq of the data sequence, in conjunction with the analytical form derived 
from our proposed one-stock model, which suggest typical features of a classical physics 
phase transition at a critical point. The double peaked form of Cq is a consequence of the 
presence of the onset of crash for the S&P500 stock index. 

Our work also differs from previous formalisms of multifractality of time series in that 
we have analytically characterized multifractal singularities and its thermodynamics inter- 
pretation. The suggested non-linear analytical forms for the supply and demand functions 
of a commodity lead us to derive theoretically the observed features of a classical phase 
transition using the simplest economic model. We believe all these novel aspects of the 
topic could stimulate further investigations on this direction and can be important to open 
beneficial discussions in the field of econophysics. 
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FIGURES 



FIG. 1. Temporal fluctuations of the Standard & Poor's 500 Index for the period 
1980-1992 as a function of trading time lags 1 < T < 220. The large, narrow fluctua- 
tions measured for T = 1 spread out on increasing T. 

FIG. 2. "Analogous" speciflc heat Cq of the S&P500 index for different time lags 
T = 1, 10, 30, 80: The + + + curve is for 1984-1988 data set, xxx for 1982-1990 and ooo for 
1980-1992. For T = 80, the full line represents theoretical results for C^+i from Eg. (P^) . 

FIG. 3. Dimension spectrum of the S&P500 index for the 1980-1992 period in the 
range g > for different time lags T = 1, 10, 30, 80, 120. 

FIG. 4. Plots of the linear p-dependence of D and Q independent of 5 (dotted lines) and 
the assumed non-linear forms for the demand and supply of Eq.(^ with \5p\ « 1 (full 
lines) . 
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